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ABSTRACT 

A pair of wave equations for the electromagnetic and gravitational perturba¬ 
tions of the charged Kerr black hole are derived. The perturbed Einstein-Maxwell 
equations in a new gauge are employed in the derivation. The wave equations 
refer to the perturbed Maxwell spinor 4> 0 and to the shear cr of a principal null 
direction of the Weyl curvature. The whole construction rests on the tripod of 
three distinct derivatives of the first curvature k of a principal null direction. 


1. Introduction 

One of the most fundamental issues of theoretical astrophysics concerns the relativistic 
sources of radiation involving black holes. The basic properties of black holes (their geom¬ 
etry, unicity and thermodynamical properties) have been discovered in the middle of the 
twentieth century. Early investigations by Petterson Petterson (1975), Chitre and Vishvesh- 
wara Chitre (1975) indicate that the minimum energy configuration of a black hole in an 
external electromagnetic field is attained when the hole has an electric charge. 

The recently discovered high-luminosity sources of X-ray radiation (Bagaroff 2001; 
Vazquez 2002) and the need for templates of gravitational radiation wave forms have am¬ 
plified interest (Hughes 2000) in models of black-hole perturbations (Wagoner 2001). It has 
been suggested (Ruffini 2001; Punsly 2001) that the observed features of gamma-ray bursters 
(GRB) can be modelled by the presence of an electrically charged black hole. There is now 
an abundant literature of this subject. Preparata and his collaborators Preparata (1998) 
present a model of GRB 971214. 

The relativistic description of radiative processes follows two distinct threads: the post- 
Newtonian expansion and black-hole perturbation theory. Despite the recent advances in 
the post-Newtonian treatment of gravitational radiation processes, it remains a laborious 
apparatus to work with. In situations involving a black hole, it would seem unreasonable to 
perturb the flat Minkowski background when the black hole metric is known to full precision. 


The first investigation of black hole perturbations has been undertaken by Chrzanowski and 
Misner (1974), Detweilcr (1977) and Chandrasekhar (1983). In the focus of their descriptions 
is the decoupled equation (Teukolsky 1973) for the perturbed curvature quantity T 0 . The 
ordinary differential equations resulting from the separation of this equation are of second 
order, with a rugged singularity structure. 

Among these pioneering studies, Chandrasekhar’s own, reviewed in Sec. 9 of his mono¬ 
graph (Chandrasekhar 1983), stands out as most detailed. In a tour de force of several 
hundred computation pages, Chandrasekhar manages to solve his seventy-six perturbation 
equations. He adds, however, the conjecture that perhaps at a later time the complexity of 
the problem will be unravelled by deeper insight. 

Chandrasekhar uses a Newman and Penrose (1962) (NP) approach both to the unper¬ 
turbed Kerr metric and to the perturbed space-time. The null tetrad vectors of the latter 
are linear combinations of those for the black-hole background. The combination coefficients 
plus the Weyl tensor perturbations are fifty real unknown functions. 

The study of electromagnetic black hole perturbations has begun with the work of 
ZerilliZerilli (1974) who considered the nonrotating system. Some years later, the problem 
of charged Kerr black hole perturbations has been taken up by Fackercll Fackerell (1982) 
and Crossman Crossman (1976). 

The purpose of the present work is precisely to help unravel the mystery of black-hole 
perturbations. 

In treatments of black-hole perturbations, it is a frequent practice to seek a suitable 
gauge fixing. The commonly applied gauge uses the Kinnersley (1969) tetrad (Chitre 1976). 
Another possible choice, promoted by Chrzanowski (1975), is the incoming (or outgoing) 
radiation gauge for the normal modes. This choice does not uniquely fix the coordinate 
gauge, however. In fact, the present work has been launched with the intent to take a 
second look at the remaining gauge freedom. We are lead to use, however a tetrad gauge 
related to the timelike Killing vector K. 

In Perjes (2002a), a simple pilot computation has revealed for stationary perturbations 
of the charged Kerr black hole that it is possible to derive a pair of wave equations for the 
electromagnetic and for the gravitational fields in a gauge related to the timelike Killing 
vector K = d/dt. We now show that this result can be extended to arbitrary perturbations 
when using the Newman et al. (1965) tetrad. To our delight, the two wave equations govern¬ 
ing the Maxwell field function $o and the shear a of a principal null direction of the Weyl 
tensor survive the generalization to arbitrary perturbations! A brief account of this latter 
result has been published in Perjes (2002b). The present paper discusses the full details. 
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In Section 2, we collect the relations describing the Kerr-Newman solution of the coupled 
Einstein-Maxwell equations. In section 3, we consider the problem of choosing the tetrad for 
the perturbed space-time, and derive a triple of expressions for three derivatives of the first 
curvature k of a principal null direction. The whole description of black-hole perturbations 
rests on this tripod. In Section 4, the integrability conditions of these derivatives will provide 
us with the key wave equations. The boundary conditions on the horizon and at null infinity 
are studied in Section 5. 


2. The unperturbed variables 

The current literature of black hole physics predominantly utilizes the Boyer-Lindquist 
coordinates ( t , r, 'd, (p ) which are related to the Kerr-Newman coordinates (t, r, 'd, p) by the 
transformation 


dt = -dt + ^^dr (1) 

dp = — dip + j^dr 

where 

A = r 2 — 2mr + a 2 + e 2 (2) 

is the horizon function. The effect of the transformation is a simultaneous reflection of the 
time and lattitude directions, and a shift of the origin of the t and ip coordinates along the 
open t trajectories and circular tp trajectories, respectively, by an amount dependent on the 
radial coordinate r. Although the Boyer-Lindquist coordinates have the salient feature that 
the two double principal null directions of the curvature are represented symmetricallyMis- 
ner (1973), they cover only the part of the original ( t,r,d,p ) coordinate domain extending 
between null infinity and the horizon. (However, this symmetry is restored in the original co¬ 
ordinates by using a twin transformation to Eqs. (1) where the directions on the trajectories 
are not reversed). 

In the original coordinates the charged Kerr metric is (Newman et al. 1965) 

ds 2 = (l - 2n ^~ e2 ) {dt - a sin 2 ddp ) 2 (3) 

+ 2 {dt — a sin 2 'Odpj {dr + a sin 2 dd/p) — (( {dd 2 + sin 2 dd<p 2 ) 
with m the mass, e the electric charge, mo the angular momentum and 


( — r — ia cos d. 
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An overbar denotes complex conjugation. The corresponding null tetrad is (Newman et al. 
1965) 


D = t 
A = n a 

5 = m a 
S = m a 


d 


dx a 

d 

dx a 

dx a 

c) 

dx a 


d_ 

dr 

1 

2 


(C + C) m - 

cc 


-1 


d_ d_ 

dr dt 


d id . . ( d d 

_ _ -U-— gj j \ iy I - — 

2 1 / 2 ( [ d'O sin d dtp \dr dt 


(4) 


The real null vector £ of this tetrad points in one of the double principal directions of the 
space-time, and the vector n lies in the intersection of the two-plane of the vector £ and 
the vector K with the null cone. (Note that the literature of black hole perturbations is 
biased towards the use of the Kinnersley tetrad Misner (1973) as a companion of the Boyer- 
Lindquist coordinates, which we do not use in this paper). In the NP notation, the following 
spin coefficients have nonzero values 


1 



1 COS '() 

2 3 / 2 ( ski'd 

m(C + C) ~ e 2 _ J_ 
2C 2 C 2C 


mb — e 2 

2C 2 C 


1 cos '0 
2 3 / 2 ( sin d 


v = ia sin d 


m C — e 2 

2‘/ 2 C 3 C 


(5) 


The nonvanishing components of the Weyl spinor are 


T 2 

^3 

T 4 


e 2 — m<C 

C 3 C 

—3iasin d 
3a 2 sin 2 d 


m C — g2 

2y 2 c“c 

m C — e2 


and the nonzero Maxwell spinor components have the form 


i 1 = 1 -F ab (rn l ’ + rffm , ’) = ^ 


$2 = F ab m a n b = 


iea sin d 

c 3 


( 6 ) 


(7) 
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In the following sections, we develop a perturbative approach to the charged Kerr black 
hole. We shall focus on the small quantities in the perturbed space-time. It is, therefore, 
important to realize that the NP quantities a, ft, e, A, 7r, r, a + /3, T 0 , Ti and <h 0 all vanish in 
the present gauge. 


3. Choice of gauge 

In this section, we introduce a new gauge for the coupled electromagnetic and gravita¬ 
tional perturbations of the charged Kerr black hole in which the wave equations decouple, 
as found in the presence of a Killing symmetry (Perjes 2002a). 

The vector £ is chosen along a principal null direction of the perturbed space-time: 

= 0 . ( 8 ) 

Unlike in the case of the principal null spinor of the exact Weyl tensor, this condition does 
not determine uniquely the spinor o A . Infinitesimal spinor rotations of the form 

o A —> o A + bi A , i A —> t A , (9) 

where b is an arbitrary first-order multiplying function, are allowed. The tetrad and spin 
coefficient transformations induced by these rotations are given in Eqs. (II) and (343) of 
Ref. Chandrasekhar (1983), respectively. 

The gauge condition ( 8 ) is preserved by the dyad transformation (9). In this gauge, the 
electromagnetic Stewart-Walker identity (Stewart 1974) 

[(D — e + e — 2p — p)(A+/i — 2y) — (5 — (3—a — 2r + 7 r)( 5 + 7 r — 2a) + crA — Kz/]<h 0 = 

— 2 [ft (A — 3q— 7 +/x) — cr (5 — 3a + /?— r) + Aft — 5 < 7 + 2 'I'i] < l>i-|-To < f > 2 - (10) 

takes a simpler form since the last term on the right vanishes. The function b is then utilised, 
following Crossman’s idea (Crossman 1976), to make the whole right-hand side vanish. The 
latter condition is employed to express the derivative 

/_ _ chu \ f A$i A 

Ak = (S - 3a + /3 - t + j a + ( 3y + 7 - n -J k - 2 Ti. ( 11 ) 

In what follows, we use a linear approximation, that is to say, we neglect products and 
higher powers of quantities which vanish for the Kerr-Newman space-time. The operators 
in the parentheses in Eq. (11) act on the first-order functions cr and ft, hence we keep their 
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form for the charged Kerr metric. For example, dq = e/ (2 1//2 £ 2 f , where the numerator 
(the electric charge e) exactly cancels in the bracketed terms containing the held dq. Thus 
we may keep our gauge condition (11) in the limiting case of an uncharged black hole, he., 
when e = 0, even though nothing remains of the original motivation for adopting it. In 
Chandrasekhar’s terms (Chandrasekhar 1983), the gauge (11) becomes a phantom gauge. 

Throughout, we express the derivatives 5p, Dp, Ap, Da, Dr, Dp and Dp from the NP 
equations (4.2k), (4.2a), (4.2q), (4.2d), (4.2c), (4.2h) and (4.2f), respectively; T>d> 2 and <5d> 2 
from Maxwell’s equations (Al) and D^>- 2 and <5T 2 from the Bianchi identities (A3) (Newman 
and Penrose 1962). 

A relation for Dk can be obtained from the NP commutator (Newman and Penrose 
1962) 

5D — D5 = (a + /3 — n) D + kA — a5 — (p + e — e) 5 (12) 

acting on dq, after expressing both Hdq and hdq from the Bianchi identities. This has the 
form 


Dk = 2 w 2 Da ~ Kp ~ 557 {D ~ ip) Dia 

+ -p— [SD + DS+(p- 4 p)S + 4 a(D + p- 2 p)] 


(13) 


A third expression for a derivative of the spin coefficient k is available from the NP equation 
(4.2b): 

5k = Da — (p + p) a + (a + 3/3 ) k. (14) 


The integrability conditions of the derivatives A k, Dk and 5k will be explored in the 
next section. 


4. Master equations 

In this section, we continue to use the gauge (11) in which Eq. (10) takes the form 

[<W — DA + (2y — p)D — 2a5 + (p + 2p)A (15) 

+ 2/ip — 2y (2 p + p) — 2 5a + T 2 + 2d>i$i] < I ) o = 0. 

Inserting the unperturbed quantities in Eq. (15), we get the first wave equation: 


Di^o — 0, 


(16) 
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where we introduce the operator 


□ s = A' s —A s+1 — + —- 

dr dr sin i) di) 


, d d , 

2 “ 1 at ~ 1 + 


s . „ a 

sm d— + s ( 1 
di) 


d 


cos 2 d 
sin 2 i) 


——2is cos d 


d_ 

dip 


+ a 2 sin 2 -2 
dt 2 


sin" i) V dp 
d 

(r 2 + a 2 ) — + (s + 2 )r + ia cos d 
dr 


d_ 

dt' 


(17) 


The spectral condition □ s d ) o = yd’o is separable in the coordinates r and i) on the basis of 
mode functions of the form f (r,i)) The homogeneous wave equation (16) can be 

brought to the form 

[(V“ + (V a + ^T a ) - 4 s 2 T 2 ] d > 0 = o (18) 

with T 2 is as given in Eq.( 6 ), s = 1 and 


F r = T^ = 0, T* = icost? + asin 2 -d cc+m(2 c c 2 ± c) ~ 2e2 , = CC- m (2C+C) +2e2 j ( 19 ) 

Equations in this class have been studied by Bini, Cherubini, Jantzen and RufhniBini (2002). 


The second wave equation arises as follows. Given the three different derivatives (11), 
(13) and (14) of k, we can derive three equations from the commutation relations. Only one 
of these equations is a second order differential equation decoupled from the other variables. 
It is obtained by applying the NP commutator 


5 A — A 5 = —vD + (r — a — (5) A — \5 — (p — 7 + 7 ) 5 (20) 

to k, and substituting the second derivatives <5A<l>i and from the NP commutators (20) 
and 

55 — 55 = (p — p) D + (~p — p) A — (a — (3) 5 — (/3 — a) 5 ( 21 ) 

for 4>i. The resulting equation has the form 


T)\(j + "D 2 $o T — 0, 


( 22 ) 


where 
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V l = 2{$ 2 $ 1 [(W -AD- 4 aS + 2 a5 + (p + p) A 
+ 7P + 4$i$i — 8aa + A p — 4 Sa + (p — 2p)p 
+ 3 (p - .7) p r- (47 — p)p + 5T 2 ] 

— 4>i [p$2 — 2z/i>i — (47 — /x)$2] -D 
+ $ 2 + 4a) - A4>! (D - 2 p)\} 

V 2 = v + 5D) - $ 2 DD] + [(p - 4 p)P + 4<& 2 $i] <M 

+ 4 [(pp — $2®i)^ 2 + ttP4*i] D + 4 [(p — 2p)P + 24>2 ®i] 
v4 = 2i> 2 4 > i [5 (3q + 7 - p) + 2(A + p + 7 - 7)7 + Pp - 4i> 2 4 > i] • (23) 


The operators T>i, X> 2 and A can be taken to have their form in the Kerr-Newman space-time 
since they act on first-order functions. The term Ak in Eq. (22) vanishes. 

It proves advantageous to introduce in (22) the ‘news potential’ ^ by 

V> = (24) 

We then have the wave equation of the form 


where 


□2^ = — J ® 0, 

Ce 


J = (e 2 — 11K,') 


d . . d id cosd\ d 

<9d dt + sin d dp sin it) dr 


- = [e 2 + m (2C - C)] 


J- r 2 
+ = e 2 

A L 


m 


( 2 C + C)] ia sind—. 


d d id cos d 

—+ msmd- + —--— 

dv dt sm v dp sm v 

d 


(25) 


(26) 


Both equations (15) and (25) pick up an additional source term in the presence of other form 
of matter, and these terms contain the corresponding stress-energy tensor T a j,. Following 
Teukolsky (1973), we expand the homogeneous solutions in quasi-normal modes with energy 
u and helicity m : 

% = f du^R (r) S™ (d) 

l,m 


(27) 
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Separating the kernel of the operator □*,, the radial function R (r) and the angular function 
S™ (d) satisfy the ordinary differential equations, respectively, 


c) d d 

A~ s —A s+1 — + 2i \(r 2 +a 2 )u—am] — + 2iu(s+2)r+2amu—A 
or o r L J O r 

1 d . d fm + scos-d\ 2 . . . 2 2 

sind——( -:— - I + (cuacosd — 1) + s — 1+A —uror 


sin i) d d d'O 


sin d 


and A is the separation constant. Solution techniques for equations of the type 
have been developed in Mano (1996, 1997). 


m (28) 

& (29) 

(28) and (29) 


5. Boundary conditions 


At the horizon, A = 0, the radial equation (28) is singular, even though the metric (3) 
remains here regular. We shall impose the boundary conditions at null infinity, r —> 00 and 
on the horizon. To this end, we perform a transformation of both the dependent variable 
and the radial coordinate: 

Z = e‘ A A’ l2 (r 2 + a 2 ) 1/2 R, = ll±£! (30) 

or A 

where Z and r' are the new radial function and coordinate, respectively. We choose the real 
function A as follows, 

A A r = (r 2 + a 2 ) u> — am. (31) 

As a result, the first derivatives of Z do not appear in the radial equation 


Z r i r i T UZ — 0 


with 


U = 


Wr [31-Tr + 2 (m — r)] — V 2 
(r 2 + a 2 ) 2 


+ [2(s + 1 )iru> + 2 amw — A —s — W] 


V — (r — m)s + i (r 2 + a~)tv — am 


W 


r 2 + a 2 


W = 


r 2 + a 2 


(32) 


(33) 

(34) 

(35) 


(i) In the neighborhood of null infinity, for large values of r, we expand the equation 
(28) in powers of 1 jr : 
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oA + ^ z „ 0 ( 36 ) 

r ) 

The asymptotic solutions of (36) and (31) are Z ~ r ±:L e ^ lu>r ' and the asymptotic form of A 
is A ~ cur. Hence the asymptotic form of the radial function is R ~ r -s e -2iwr' an( ] r ~ s ~ 2 . 



(ii) At the event horizon, r = r + or r' —> —oo, the radial equation becomes 

V 2 


Zyr ' (a^ + rl) 2 


Z Rj 0. 


(37) 


The boundary conditions on the event horizon will be selected by taking into account 
the absence or presence of incoming or outgoing radiation. 

The boundary condition for the angular equation (29), imposed at i) — 0 and at d = n 
is that S must be regular on the boundary. As before (Teukolsky 1973), this leads to the 
Sturm-Liouville eigenvalue problem and yields a complete set of complex solutions. 


6. Ramifications 

The Maxwell equations for electromagnetic perturbations on a fixed space-time are 
homogeneous. The circularly polarized normal modes have the structure (j) A (r, $) e l ( m ‘P- u ’ t )_ 
The interaction with the gravitational field, via the stress-energy tensor 0 A 0 B , introduces 
a mixing among the electromagnetic modes, even in the linearized theory. Thus, the full 
electrovacuum perturbation equations cease to be homogeneous in the modes. In quantum 
field theory, perturbative methods are available for such interacting systems. Strangely, no 
comparable treatment for the corresponding classical system is known. This regrettable 
backlog in the classical theory, however, will help illuminate the origin of the long-standing 
difficulties in finding a separable equation for the perturbations. 

In this paper, the following picture emerges for the classical electrovacuum pertur¬ 
bations. There exists a subset of the field equations, consisting, in part, of gravitational 
equations and of some of the Maxwell equations not containing any mode mixing. (This is 
because the complex conjugate electromagnetic stresses are absent from these equations). 
From these relations alone, it is possible to obtain a pair of wave equations for the quantities 
$o an d er. Thus a normal mode expansion for this doublet of fields is available. 

Given the perturbed quantities <f>o and cr, Eq. (13) is an inhomogeneous first-order ordi¬ 
nary differential equation for the function k. Integration of Eq. (13) yields the r dependence 
of the spin coefficient k. Once this is known, the r dependence of expressions like An is 



fixed. We then get 'hi by solving the linear algebraic Eq. (11) for \&i. Mode mixing occurs 
only in the further gravitational and electromagnetic perturbation components. 
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